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ABSTRACT 

In the case of gravitationally-lensed quasars, it is well-known that there is a 
time delay between occurrence of the intrinsic variabilities in each split image. 
Generally, the source of variabilities has a finite size, and there are time delays 
even in one image. If the origin of variabilities is widely distributed, say over 
£ 100 pc as whole, variabilities between split images will not show a good 
correlation even though their origin is identical. Using this fact, we are able to 
limit the whole source size of variabilities in a quasar below the limit of direct 
resolution by today's observational instruments. 

Subject headings: galaxies: active — galaxies: nuclei — gravitational lensing - 
quasars: individuals (Q0957+561) 
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1. INTRODUCTION 

Since Liebes (1964) and Refsdal (1964) have reported meaningful aspects of gravitational 
lensing phenomenon, many researchers rushed into the field of gravitational-lensing study, 
and presented many interesting results. This situation is not altered in these days. 

One of the most interesting gravitational-lens phenomena is quasar lensing. This is 
caused by a lensing galaxy (or galaxies) intervening observer and quasar. In the context 
of cosmology, it will be possible to estimate Hubble's constant from a time delay of the 
quasar variations between gravitationally-lensed, split images. The most successful study is 
by Kundic et al. (1997, hereafter K97). They monitored Q0957+561 for a long time and 
performed robust determination of the time delay. From their own result, they evaluate 
Hubble's (H ) constant as 64t^ km s _1 Mpc -1 based on the lens model constructed by 
Grogin and Narayan (1996, hereafter GN). 

On the other hands, concerning the structure of quasar, we will discriminate the 
structure of central engine according to the effect of a finite source size. Recently, Yonehara 
et al. (1998, 1999) performed realistic simulations of quasar microlensing, and showed that 
mult i- wavelength observations will reveal the structure of accretion disk believed to be 
situated in the center of quasars. Furthermore, using precise astrometric technique, Lewis 
and Ibata (1998) indicated that it is also possible to probe the structure of quasar from 
image-centroid shift caused by microlensing. Observationally, in the case of Q2237+0305, 
Mediavilla et al. (1998) detected a difference between an extent of the continuum source 
and that of the emission-line source by two-dimensional spectroscopy, and limit the size 
of these regions. Thus, quasar-lensing phenomena are a useful tool to probe not only for 
cosmology but also for the structure of quasar. 

Following these interesting researches, we propose a method to estimate, in this letter, 
the effect of a finite source size on time delays of the observed quasar variations between 



-4- 



each gravitationally-lensed, split image, and to judge whether it is negligibly small or not 
and to limit the whole size of the source of quasar variability. This is important because 
no such limitation has been done yet although the size of each variation, "one shot", 
had already been obtained order of days assuming causality in the individual source of 
variations. 

In section 2, I describe the basic concept of this work, and simply estimate the time 
delay difference. Next, I present some results of calculation for the case of Q0957+561 in 
section 3. Finally, section 4 is devoted to discussion. 

2. BASIC CONCEPT 

The basic idea that we wish to present in this letter is schematically illustrated in 
figure [I]. Suppose the situation that a quasar is macrolensed by lensing objects so that its 
image is split into two (or more) images. The angular separation between these images 
is large enough to observe individually, say apparent angular separation is £ 1 arcsec. 
If we observe such quasar images, we will realize the intrinsic variabilities of quasar in 
each image as in the case of an ordinary, not gravitationally lensed quasar (e.g., recent 
optical monitoring results are shown in Sirola et al. 1998). Because of the macrolensing 
effect, generally, the variabilities in such a quasar are not observed in both images at the 
same time. There is a time delay between these quasar images caused by a light path 
difference from the light path without lensing objects which originates from gravitational 
lens effect (e.g., see Schneider, Ehlers, & Falco 1992, hereafter SEF). These facts are nicely 
demonstrated by K97. 

However, previous studies related to the time delay caused by gravitational lensing 
were not so much concerned with the source of variabilities, and the source of variabilities 
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was treated as a point source. This treatment is reasonable, if the whole source size is 
negligibly small compared with the typical scale length over which a time delay changes. In 
contrast, actually, we only know that the source of quasar variabilities is smaller than the 
limit of the observational resolution, say £ 1 arcsec (e.g., in the case of HST observations, 
~ 0.1 arcsec), and we do not know whether the whole source size is small or large compared 
with the scale length over which a time delay changes. Therefore, first, we should try to 
consider the effect of a finite source size on the expected observed light curve in quasar 
images. 

Then, if we include such an effect, what do we expect to see ? The answer is easily 
understood by figure [T]. For simplicity, I consider only two images (image A and B) of 
the lensed quasar, and the source exhibit only two bursts ("burst 1" and "burst 2", they 
occur in this order on the source plane) with some time interval (Atburst)- The origin 
and separation of such bursts are not specified, we assume that these two bursts are not 
physically correlated, in other words they appear randomly. Additionally, we set a time 
delay difference between the position of the "burst 1" and the "burst 2" on image A as At a 
and that on image B as At-g. 

In the case of Atburst ^ |A£a — Aiel, light curves of two images show apparently very 
similar feature, instead of its time delay at the very center. Although the shape of light 
curves is altered from intrinsic one by the effect of finite source size as is depicted in lower 
left part of figure [I], we can easily identify these two light curves are intrinsically the same 
one. Thus, we are able to obtain a robust time delay between two images. 

On the other hands, in the case of Atburst ~ |A£a — Aisl, a previous fact does not hold 
any more. In this case, time interval between two bursts is significantly modified by the 
effect of its apparently large time delay difference (|A£a — Aie|)- In such a situation, we 
can no longer conclude that light curves from two images have the same origin, even if we 
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include an effect of time delay for the case of point source. We may seek for the reason for 
this to microlensing or something exotic. In other words, there will be no good correlation 
between light curves of two images. This is a serious problem not only to determine time 
delay or Hq but also to construct a quasar structure, to determine the origin of variabilities, 
or some other problems. 

Here, I will make a simple estimate of time delay difference between different parts of 
the source, i.e., the effect of a finite source. In this estimate, I define (3,6a, 6b as angular 
positions of the source center and those of the centers of two images. Therefore, 6a and 6-q 
are the solutions of well-known lens equation (e.g., SEF), 

(3 = 6-a(6), (1) 

where, a is a bending angle caused by intervening lens object (s), i.e., gravitational lens 
effect. Furthermore, time delay from un-lensed light path (r) in the case of the image 
position is 6 and the source position is (3 is written as 

t = { -±±^V \6-(3\ 2 - Q±£*L*(e). (2) 

2c c d 

Here, z Q \ is redshift from observer to lens, T> is effective lens distance that by using angular 
diameter distance from observer to lens {D \), from observer to source (D os ), and from lens 
to source (As), written as V = AiAs/ As, and \& is so-called "effective lens potential" 
(e.g., SEF). Insert each image position into equation (0) and subtract one equation from 
the other, we obtain the well-known time delay expression between images A and B (Atab), 

A tab (P) = {\0a-P\ 2 - \9 b - (3\ 2 ) - ii±^{¥(6l A ) - (3) 

Additionally, if we assume the position that is offset by d(3 from the center of the source 
and write d&A and dd-Q as image positions from the center of the image, these variables 
should fulfill the lens equation (|l|) again, i.e., 

(J3 + d/3) = (6i + d6i) -a(6i + ddi), (i = A or B) (4) 
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or, subtracting this from equation (|l]) and adopt Taylor expansion to a, we obtain another 
expression of equation (f|), d(3 = d6i — S/ga(9i)d9i + ■ • •. 

Subtracting At ab (/3) from Ar AB (j3 + d/3), c.f., equation (|3|), and using equation (|l|) and 
equation (||), we are able to obtain the time delay difference between the center of source and 
the other position offset by d/3 from the center of the source (Sr AB = Ar AB ([3+d{3) — At ab ({3)) 
on the source plane. 

Moreover, by definition of effective lens potential, bending angle a is related to the 9 
through the derivative of effective lens potential ^(9) as V 'e^(9) /T>c 2 = a(9). Since we 
are considering the origin of quasar variabilities, the source size is ~ kpc at most and the 
distance from observer is typical cosmological scale ~ 1 Gpc. Thus, its apparent angular 
size is ~ kpc/lGpc = 10~ 6 rad. This seems to be small compared with image separation 
and the scale of bending angles which is typically a few arcsec. For <5tab, accordingly, 
we can adopt a Taylor expansion to a and neglect the higher terms than first order 
assuming d[3 -C (3 and d9i <C 9^. After using some algebra and putting R as actual 
off-centered distance on the source plane, i.e., R = \dj3\ ■ D os , we are able to evaluate time 
delay difference as follows, 

5t ab ~ ^ + 2 ^ V{2a(9 A )V e a(9 A )d9 A - 2a(9 B )V e a(9 B )d9 B } 

_ (l+£ol) _ V ^ b)c ^ b | 

cr 

(1 + z ol ) 



-V [{V,a(fl A ) - l}a(9 A )d9 A - {V e a(9 B ) - l}a(9 B )d9 B ] 

c 



c 



-V\(9 B -9 A )-df3\ (5) 



This one-dimensional evaluation is somewhat overestimated, however, for the calculations 
above, we did not use any restriction about lens model, and equation (||) seems to be 
appropriate for any lens models and lensed systems except in some special situations, e.g., 
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in the vicinity of caustics (or critical curves). 

Consequently, considering the fact that quasar optical intrinsic variabilities have 
timescale day ~ month, equation (|6]) indicates that correlation between light curves of two 
images shown, in worst cases, will disappear, if the origin of quasar variabilities is extended 
over ~ 1 kpc, i.e., maximum off-centered burst occurs at ~ 1 kpc from the center of quasar. 

3. EXAMPLES OF Q0957+561 

Finally, we will show some impressive result for the case of Q0957+561 which is the 
first detected lensed quasar by Walsh, Carswell, & Weymann, (1979). 

To demonstrate how the extended source effect works on the time delay determination 
in an actual lensed quasar, here, I will present simulation results of Q0957+561 as one 
example. Using equation (^|), we are able to estimate a time delay difference between 
same source positions at different lensed images. In this CclSG, clS IS well known, if we use 
z Q \ ~ 0.36, z os ~ 1.41, \6 B — 9 A \ ~ 6" (e.g., GN) and assumed that H ~ 60km s _1 Mpc -1 , 
we will obtain Stab ~ 50 day for the source with a size of 1 kpc ! 

Furthermore, to obtain more realistic results, we used isothermal SPLS galaxy with 
compact core as an example of lens model for Q0957+561 (details are shown in GN), 
adopted parameters listed in table 7 in GN as "isothermal" SPLS and calculated time 
delay difference between images center and off-centered part of images (Stab)- For 
this calculation, we set Q = 1.0 for simplicity and took on convergence k ~ 0.22 and 
H ~ 64 km s _1 Mpc -1 to reproduce the observed time delay following K97. The resultant 
time delay contours compared with the image centers on the source plane are depicted 
in figure |2]. On image A (left panel), a gradient of the contour is almost in the negative 
y-direction, although that of image B (right panel, this time delay advanced ~ 417day ) is 
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almost in the positive y-direction. Additionally, time difference between the same position 
on the source reaches order of ~ months for the case of the source with a ~ kpc size, 
therefore, we expect disappearance of correlations between the light curves of image A and 
that of image B. Here, from equation (|5|), we can easily understand why do the contour 
lines show almost straight and perpendicular to y-direction. Product of (6* a — 6b) and d/3 
in equation means that time delay difference determined mostly by the element which 
parallel to (8 a ~ #b) of displacement d/3. Therefore, the time delay difference significantly 
alters along the (6 a — #b) direction and almost constant along the perpendicular direction. 

Moreover, we simulated expected light curves of variabilities in both quasar images 
using superposition of simple bursts with triangled shape and duration of ~ 10 days which 
are randomly distributed in time, in space, and in amplitude. For the whole source size, 
I consider three cases, 1 kpc, 100 pc and 10 pc. Using the same procedure to produce 
figure 0, we calculated time delay from the center of image A over both images, randomly 
produce bursts, sum up all bursts and finally obtain expected light curves as presented in 
figure "Residual light curves" produced by subtracting properly-shifted light curve of 
image B from that of image A, are also shown in the figure. In the case of the smallest 
source, R = 10 pc, and still in the case of middle source size, R = 100 pc, we can easily 
recognize the coherent pattern in the light curves of images A and B but with time delay of 
~ 417 days advanced light curve of image B. Time delay between two image centers is able 
to be determined fairly well. However, in the case of largest source, R = 1 kpc, it seems 
no correlation between two light curves even if we already know the time delay between 
two image centers, and we may misunderstand that the variabilities did not originate in 
source itself ! This feature is far from the observed properties that the time delay between 
two images is determined easily even if we fit them by eyes. Therefore, I conclude that the 
size of source that is origin of quasar variabilities should be smaller than 100 pc, namely, 
maximum acceptable size is order of ~ 10 pc from this simple simulation. 



4. DISCUSSIONS AND COMMENTS 



As we examined, if we include the finite source-size effect to the time delay 
determination from quasar variabilities, correlation between expected light curves of each 
lensed image will disappear in the case of the size is sufficiently large, say ~ lkpc. Using 
this fact, we can limit the size of the region where quasar variabilities are produced, 
from the correlation between light curves of multiple lensed quasar images. Furthermore, 
since the size of the origin of intrinsic quasar variability reflects a physical origin of the 
variabilities, we can also determine the origin of the variabilities, e.g., whether it is disk 
instability (Kawaguchi et al. 1998) or star burst (Aretxaga, Cid Fernandes, & Terlevich, 
1996). Particularly, in the case of Q0957+561, the origin of variabilities has a size smaller 
than 100 pc. This value is consistent with disk instability model, because of its small size 
(~ O.Olpc for 1000 Schwarzschild radius accretion disk surrounding 1O 8 M supermassive 
black hole). Starburst model can be rejected, since starburst region is 100 pc ~ 1 kpc. 
Hence, for the origin of intrinsic quasar variabilities, the disk instability model is more 
preferable, as was indicated by Kawaguchi et al. (1998) already. To draw this conclusion 
more critically, we should do this study more precisely in future. 

Additionally, the fact that a larger source size tends to reduce a good correlation 
between the light curves of each image provides an answer to the question why time delay 
determination from radio flux gave a wrong answer except recent works, e.g., Haarsma et al 
(1999). Generally, radio emitting region is believed to have a larger size than that of optical 
photon because of the existence of large radio lobe and/or jet component, and the effect we 
have shown in this letter may be significant. Thus, robust determination of the time delay 
seems to be difficult. 

If such a effect is significant in the well-known lensed quasar Q2237+0305, microlens 
interpretation of individual variabilities (e.g., see Irwin et al. 1989) will be rejected. 
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Fortunately, however, this may not be the case because for this source, caused by its quite 
nice symmetry of lensed image, the effect seems not to be so significant and intrinsic 
variabilities will be expected to appear in every images with good correlations. 

If we develop this technique furthermore, and adapted to another multiply-imaged 
lensed quasar, we will determine the size of the most interesting part in quasar. 

The author would like to express his thanks to Toshihiro Kawaguchi, Jun Fukue for 
their number of comments, and to referee for his/her meaningful suggestions. The author 
also grateful to Shin Mineshige for his valuable discussions and kind reading of the previous 
draft. This work was supported in part by the Japan Society for the Promotion of Science 
(9852). 
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Fig. 1.— Schematic view of the finite-source size effect for time delay between 
gravitationally-lensed, split quasar image (see section 2). 
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Fig. 2. — Contour of time delay difference at image A (left panel) and image B (right panel) 
compared with their own center. The y-axis positive to north same as GN and the axis is 
almost parallel to the 9a — #b direction. Source center is shown by filled circle at the center 
of each figures and long-dashed circles present a circle with a diameter of 1 kpc, 300 pc, and 
100 pc from outer one to inner one, respectively. The time delay contours of leading parts 
(solid lines), preceding parts (dotted lines), and no time delay from the center parts (dashed 
lines) are presented. Concrete values are also shown in the figures. 
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Fig. 3. — Simulated light curves (upper one of each panel) of image A (solid curve) and 
B (dashed curve), and its residual light curves (lower one) for some different source size. 
Numeric values on the horizontal line are fluxes normalized by a standard deviation of light 
curve of image A (<7a). The whole source sizes are 10 pc (upper panel), 100 pc (middle 
panel), and 1 kpc (lower panel). The curves are arbitrary shifted in the vertical direction. A 
light curve in dotted box of image A and B should be the same in the case of point-source 
treatment. 



